High-momentum asymtotics from the 
-Q^-- Fock-Feynman-Schwinger path integral 



Yu.A.Simonov 



5-H 

CL|i Institute of Theoretical and Experimental Physics 

■ 117218, Moscow, B.Cheremushkinskaya 25, 

^ ■ Russia 



Abstract 

cn 

'sj" I The asymptotics of n-point Green's function at large external mo- 

menta is obtained in the exponentiated form using the Fock-Feynman- 
, Schwinger representation for propagators in the external field. The 

I method is applied to gauge theories such as QCD and QED, and the 

• Sudakov form-factor is calculated as an example in QED and me- 

son form-factor in QCD. Nonperturbative contributions can be conve- 
Qh! niently included, as it is demonstrated in the example of the confine- 

^ ' ment correction to the form-factor. 

> 

X 1 Introduction 



The high-momentum asymptotics of QCD and QED amplitudes is under ac- 
tive study for several decades [1,2]. The basic method used is the summing 
up of the series of dominant Feynman amplitudes keeping the leading asymp- 
totic term in each of them. These terms appear to be of the type a^ln^Q'^ 
and in case r = 2n one has the so-called double logarithmic asymptotics, 
which e.g. defines behaviour of the QED form-factor at large [3]. An 
extension of the double-logarithmic asymptotics to other processes has been 
made both in QED [4] and QCD [2,5]. 

One of attractive features of these results is the exponentiation of leading 
logarithmic terms coming from each Feynman amplitude due to summation 
over n, which might be a hint, that there can exist an alternative, more direct 
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method of deriving this result. A step in this direction is made in the present 
paper. We start with a most general n-point Green's function, depending on 
n external momenta pk-,k = l,...n. Two kinds of theories are considered: 
QED and QCD, but the generalization to the EW case is straightforward. In 
this first study we confine ourselves to the simplest setting of the problem: 
all combinations of momenta are large compared to masses, external Mj and 
internal mj, 

pI->MI ml I, k=l,...n (1) 
Moreover the following condition is assumed to hold 

\PiPk\^Pi,Phhk^l,...n (2) 

Our main object of the study will be the one-fermion-loop amplitudes which 
assumes the limit of large Nc for QCD. However one may expect that the 
inclusion of additional fermion loops gives only small correction to the leading 
asymptotics, so that the one- fermion-loop amplitudes are dominant in the 
asymtotics (1) and (2) both for QCD and QED. 

For QCD there is a special problem of nonperturbative (NP) contribu- 
tions, which can be handled as in [6], i.e. dividing the total vector field 
Afj_ = + a^, where is the NP background and is the perturbative 
field. However with the conditions (1) and (2) the standard expectation is 
that the NP contribution dies out as a power, e.g. cx/p^, i;^p^y where a is 
the string tension (or any equivalent combination like gluonic condensate) 
made of the NP field i?^. We shall argue in the last section of the paper that 
actually NP effects can be also exponentiated and may strongly modify the 
perturbative result. The plan of the paper is as follows. In section 2 we intro- 
duce FFSR for the n-point Green's function . In section 3 the approximate 
path integration is done in a way pertinent to the high momentum kinemat- 
ics (1),(2) and double logarithmic terms in the asymptotics are calculated. 
In section 4 the Sudakov vertex asymptotics is identified and formfactors 
are discussed. In section 5 the influence of confinement is considered and 
as a result the confinement correction is calculated in addition to the DLA 
term. Discussion of other applications is given in the conclusions, and in the 
appendix the spin-dependent terms (crF) in the fermion propagators are es- 
timated and shown to be subleading. Note, that the situation with NP fields 
changes, if one abandons condition (1), since the Green's function under con- 
sideration has bound state poles at pj = Mf, and those are mostly due to the 
NP forces between quark and antiquark. As it is, however, both conditions 
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ensure that to the leading approximation one can consider only perturbative 
(gluon or photon) exchanges between fermion Unas and the main point is 
how to sum up those in the most effective way. 

To this end we shall exploit the Fock-Feynman-Schwinger representation 
(FFSR) for the fermion propagator in the external field [7,8], which was used 
previously for similar purposes in [9]. The path-integral method have been 
used for the calculation of the Green's function asymptotics in QED for a 
long time, see e.g. [10] and refs. therein. More recently the so-called world 
line method was introduced in [11]. Below we follow the formalism of [9], 
which allows for a more economical treatment of spin. 



2 Integral representation for the fermion prop- 
agator 

For the spin-1/2 particle propagator in the external vector field one has 
the FFSR (we start with the Euchdean metric for convenience) 

S{x, y) = {m-D) ds{Dz),ye-''^^{x, y) (3) 

J 

where notations used are 

K^mh+-J^ dri^r, D,^d,- igA„ (4) 

(5) 

PX PS 

^a{x,y) ^ PAPpexpig A^dz^j,expg dra^^^F^,^, (6) 

Jy Jo 



and 



t^Mf = ^(7/.7i^ - lulfi), ^{n) = z{n) - z{n - 1). 



Consider now the n-point Green's function with external momenta at n ver- 
tices equal to pi (see Fig.l) 

G{pi, ...pn) =< J\{Pi)-JniPn) >, Ji{x) = '4^{x)Tii^{x) (7) 
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Insertion of (3) into (7) for the one-fermion loop yields 



^=l 

(8) 

Wc shall disregard in what follows the factors Ti{mi — Di) since we shall be 
interested only in the exponentiated contributions; doing the dx^'^^ integrals, 
one obtains 

n 

1=1 

where 

<5» = / n n 15'""' < r^f"" - > 

and 

i=l 

The integral d'^q denotes the integral over one of qi, all others being expressed 
through it and all pi. 

We note that < W^- > is a gauge invariant quantity summing all the 
perturbative exchanges inside the fixed Wilson contour, defined by the set 
{^^^\k)}. In addition to the usual Wilson (charge) vertices, there are also 
magnetic moment vertices aF, hence the notation < W^j >. 

We concentrate now on the contribution of the Aj^ field in (6), referring 
the reader for the discussion of the aF term to the Appendix. 

The crucial step for what follows is the use of the cluster expansion the- 
orem [12], which yields for < Wa >^<W> 

°° (ioY 

r=l ki 

(10) 

Here double brackets denote cumulants [12], the lowest order contribution 
(in the exponent) is expressed through photon (gluon) propagator, which in 
the Feynman gauge is (the gauge is irrelevant since < W > is gauge invariant) 

< A,(.)AAz') >= (11) 
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Here C2(/) is the quadratic Casimir operator for fundamental representation, 
1 is the unit color matrix, for QED one should replace C2I — * 1. 

In what follows we confine ourselves to the lowest contribution (11) in 
(10) and show that it yields the double logarithmic asymptotics, leaving 
next order terms for the future. 

First of all one can persuade oneself that the approximation (11) yields 
in (9) all diagrams with exchanges of photon/gluon lines between fermion 
lines, all orderings of lines included. For QCD this means the following: 
all orderings, i.e. all intersection of gluon lines in space-time are included, 
except that color ordering of operators is kept fixed. Since commutator 
of any two generators t"' is subleading at large Nc, it means that (11) sums 
up all exchanges including intersections of gluonic lines in the leading TVc 
approximation (c/ the discussion of this point in [13], note that e.g. the 
nonplanar diagram 0{g^) is given by the quartic cumulant in (10) and the 
latter is suppressed by the factor 1/Nc). 



3 The Gaussian integration 

Our next point is the integration over d^{k) in (9) which is Gaussian in the 
main term Ki, defining the measure of integration; therefore we shall do 
it expanding the exponent in (9) around the stable fixed point C,, which is 
obtained by differentiating the exponent in (9) with respect to ^^''\k). One 
has 



47r2 j^,{z(^{k)-z:>{k'))^ 

Here e.g. z^'^^k) = Z]j=i Y.t=i ^'^^^i^), where we have chosen as the origin the 
coordinate x^^^ of the first vertex, and all other coordinates are calculated 
using the connection x^'^ — x^'^~^^ — J2k=i^^^K^) with the cyclic condition 

r^{n+l) _ ^(1) 

One can solve the system of equations (12) iteratively expanding in powers 
oi g^, the first two terms are given in (12), where one should replace ^^'^^ inside 
the curly brackets by 2ieiq^'^\ If one represents the exponential appearing in 
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(9) after insertion of (10) as exp{—f{^, q)), then one can write 

(13) 

The Gaussian integration in (9) finally yields 

G'n- /T^n^^.e-/^^""')-^*'-^"- (14) 
l^'^J i=i 

where the matrix (p is 



i2 



i- 1 d' 
'^'^'^ " 2ae«(A;)ae(^)(n)^^^' '^^ 



(15) 



The most important for what follows is the term f{^,q) which can be 
written as (at this point we reestablish Minkowskian metric) 

where we have defined — ksi, and 

A^,=i:w\i<i (17) 

The integral in the last term on the r.h.s. of (16) can be written as 
SiSj(gWg(i))/(s^ g)^ where 

/■I /■! dadf] 
UAs,q) = I I (as,g»-/?s.g(i)-A,,)2 ^^^^ 

The diagonal terms, lij, with i = do not contribute to the asymptotics and 
contain only selfenergy divergencies, which are of no interest to us in what 
follows. Therefore we shall consider nondiagonal terms with i ^ j. 

Let us first study the term with i = j — 1 ('the dressed vertex contri- 
bution"), and Aj_j+i = see Fig. 2. Then (18) is reduced to the form 
which will be studied below 

1 /-i dadj /-i /-i d(xd:3 



{asiq^) + (3sjq^i)Y Jo Jo {a^a^ + (3%'^ + 2al3{ah) 

(19) 
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with a = Siq^'^\ b — Sjq^\j = i + 1. As it stands the integral (19) diverges at 
small a, (3 (or at small Tj.Tj in (16)). The origin of this divergence becomes 
physically clear, when one expresses the distance z^'^^ from the vertex position, 
(we go over to the Minkowskian space-time) 

^W^2g«Ti,2(^) = 2g(^V^-. (20) 

The quasiclassical motion (20) cannot be true for small when quantum 
fluctuations wash out the straight-line trajectories, and the lower limit Tmin 
can be obtained from the quantum uncertainty principle 

A^Ag ~ (^« - z^^^){q^^ - g(^)) ~ 1 (21) 

We shall be interested in the kinematical region, where 

|(^W^W|»(^W)2^(^0-))2^ (22) 

and Tmin then is found from (21) to be 

1 , , 

^""''^ ^ 2p^ ^^^^ 

Using (23) one can easily calculate the integral (19) since the term 2aP{ab) 
in the denominator of the integrand in (19) always dominates. The result is 

= 2,^,^^^(g(.)g(.+i)) ^^(2(gg0^i)^^(2(ggO^^+i) (24) 

The integration of the general term lij with j ^ i — + l can be done using 
the expressions for the Spence functions. However in the general case the 
lower limit Tmin is inessential and the double logarithmic situation does not 
occur unless there is a large ratio, \j3i3hl\ ^ i We leave the detailed study 
of this point to the future. 

4 The 3— point Green's functions 

We start with the open triangle. Fig. 2 corresponding to the Sudakov vertex 
function asymptotics, 

Gs = {-iq + m)-^r{q, q'){-iq' + m)"^ (25) 
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In this case there is no integration over dq in (9) and only one integral J12 
(18), (24) is present (we disregard as before the selfinteracting pieces which 
do not contribute to the asymptotics). 

Insertion of the integral (24) into (16) and integration over dsids2 in (14) 
which yields in (24) (with logarithmic accuracy) replacement — > ^, finally 
leads one to the answer for QED (C2 = 1) 

VM) ~ eM—ln'Mlln^) (26) 

which coincides with the known Sudakov asymptotics [3]. 

We turn now to the case of QCD, where the basic triangle diagram is 
closed due to color gauge invariance and try to find out whether the kine- 
matical region (22) plays important role in the integral over dq in (14). 

In the general case, when all qi are unconstrained and expressed through 
three external momenta Pi,P2,P3 and one integration variable, the region 
(22), yielding double logarithmic asymptotics (DLA) (24), is suppressed due 
to large values of f{C,q) in the exponent. As a result the integral over dq 
does not lead to the DLA form for G3. 

The situation changes however, if one considers instead of G3 the form- 
factor, i.e. when the pole terms are factored out from the vertices 2 and 3 
see Fig. 3, and vertex functions appear there. To simplify matter, one can 
consider for the formfactor the same representation (14), where under the 
integral one introduces vertex functions ipi{ki),i — 1, 3, where 



The definition (29) yields in the c.m.system of particle 1 or 3 the familiar 
relative momentum of two emitted fermions. The presence of ipi imposes 
restriction on momenta qi, namely 

klkl<K' (28) 

where /t^ is some hadronic scale. 



We can define in the Breit system momenta as follows p^^\ p^'^\ p^^'^ 
-f ), (0, Q), {Po, f) where > and p^ ^ + f . 
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One can easily see, that (30) constrains the region of integration over 
dq = d^qi to the region |g| ~ |?o ~ Pol ~ i^i a-nd conditions (22) are 
satisfied. 

Hence in this case one can reproduce the Sudakov asymptotics (26), where 
Igg'l ^ ^^,q^r^ g'^ r^K\ a,C^ (29) 

5 Modifications due to confinement 

As was discussed in the Introduction, in the kinematical region (1),(2) one 
usually assumes that all distances in the n-point function are small and 
therefore nonperturbativc interaction is ineffective. Here we show that this 
is in general not true, and in particular for the DLA processes large distances 
appear naturally. Indeed, the distances between two straight-line trajectories 
(20), starting from some vertex, I = \z — z'\ ~ 2\qT — q'r'l may become 

large, I = 0{^^jj), at the end points, r = s ~ ^, r' = s' ~ ^. The 
nonperturbative background produces the area-law term for < W^. > in 
(9), < Wcr >~ ea;p(— (xS'min)) where .Smin is the area of the minimal surface 
bounded by the line of trajectories, {^''''\k)}, and the natural measure for 
infiuence of confinement is the value of crSmin foi' trajectories (20). For the 
triangle formed by the latter the area is Smin ~ \J {zz'Y ~ |??'|tt'. It 
is clear that confinement starts to play role at the value r ~ r' ~ tq = 
((T|gg'|)~^/^ and at later proper time trajectories are no more straight lines. 
To estimate corrections to the trajectories one can add the term crSmin to the 
effective action (13). Differentiation of aSmin with respect to ^j^^\k) yields an 
additional term on the r.h.s. of (12) ~ ar^, where r-^ is the component 
of = 2i{q^T — q^f) perpendicular to Since ti < Sj, Si ~ ^, one 

obtains correction O(^) to the leading term q^^^ in the curly brackets of (12). 
Hence the condition of stability of Sudakov asymptotics with respect to the 
nonperturbative (confining) effects is qf ^ a. 

However there is nevertheless an additional contribution due to crSmin in 
the exponent, which obtains when one inserts in Smin straight-line trajecto- 
ries (20). With the conditions (22) one gets 

< ex.p{-aSmin >eik^ ea;p(-^^^^) (30) 

I I 
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One can notice that the term in the exponent is not necessarily small and 
may substantially modify the double-logarithmic asymptotics. 

6 Conclusion 

The formalism described above is a direct development of the relativistic path 
integral, using the proper time, which we call Fock-Feynman-Schwinger Rep- 
resentation (FFSR) [7,8]. As in other versions of the path integral method 
[10,11] it is convenient to yield results already in the exponentiated form. 
We kept here the bispinor variables in the original Feynman form (the term 
< Wa > in (9)) which in principle allows to calculate spin-dependent ef- 
fects also in an exponentiated form, however it was shown in the appendix 
that those effects are subdominant. One of the important advantages of the 
present formalism is that it allows to take into account NP effects also in the 
exponent, as it was demonstrated in section 5. 

There are possible many lines of development, including application to 
Drell-Yan and heavy-quark production, as well as finding connections with 
the factorization theorem, as it was done in a similar formalism in [14]. 

The author is grateful to A.B.Kaidalov for useful suggestions and to 
L.N.Lipatov for discussion of the results. 

7 Appendix. Contribution of the magnetic 
moment term a^uF^^y. 

We prove below that (cxF) terms do not contribute to the DLA, but yield 
the subleading terms. 

The terms to be estimated are 



To illustrate the procedure it is enough to consider the O(gf^) terms, namely 



where Dp^^o-i/ =< FpiJ.{z)F^v{z') > can be written perturbatively as {u = 




(A.l) 




(A.2) 
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Id 1 
Dp^,av = ^[^K'^/''^ - '^^^1^) +P^^ (A.3) 



Commuting 7~matrices (A. 2) can be rewritten as 

K47A.-7A^)/ dr dr'^ (A.4) 
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In the lowest saddle-point approximation for ^, Eq. (12), one can replace 
according to (20) 

2i{q,T + q'T') (A.5) 



At the same time the lower limit of integration over r, r' in (A.4) should be 
replaced by Tmin — ijgTi' integral can be estimated as 

' T^^h + (?,gt + i?.)^2 + ^/3] (A.6) 



2\qq 

where /„, n = 1, 2, 3 grow at most logarithmically at large \qq'\ and all terms 

2 I'Z 

of the order Ofr^, -r—n) have been neglected. Thus the terms aF do not 

V \qq'\ ' Igg'l^ 

contribute to DLA in the order O(gf^). One can follow the same procedure 
to higher orders and persuade oneself that leading powers of logarithms, 
aln'^\qq'\, do not appear. 
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